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Abstract 

The one-way measurement model is a framework for universal quantum computation, in which algorithms are partially 
described by a graph G of entanglement relations on a collection of qubits. A sufficient condition for an algorithm to perform a 
unitary embedding between two Hilbert spaces is for the graph G, together with input/output vertices I,0 V{G), to have a 
flow in the sense introduced by Danos and Kashefi For the special case of |/| — \0\, using a graph-theoretic characterization, 
I show that such flows are unique when they exist. This leads to an efflcient algorithm for finding flows, by a reduction to 
solved problems in graph theory. 



I. INTRODUCTION 

The one-way measurement model is a framework for 
quantum computation, first presented in P, 01 ■ Trans- 
formations of quantum states in the one-way measure- 
ment model are essentially described by a sequence of 
single-qubit measurements (where the choice of measure- 
ment may depend on earlier measurement results in a 
straightforward way) performed on a many-qubit entan- 
gled state. The many-qubit system includes some num- 
ber of input qubits I in an unknown initial state, and a 
collection of auxiliary qubits prepared in the |-|-) state. 
These are operated on by a collection of entangling oper- 
ations, described by a graph G whose edges uv G E{G) 
are pairs of qubits operated on by two-qubit controUed- 
Z operations. After the sequence of measurements, any 
qubits left unmeasured still support a quantum state, and 
are interpreted as an output system O. 

Algorithms in the one-way model may be obtained by 
translating from the circuit model: we may decompose 
a unitary operation U into one- and two-qubit unitaries 
which have known implementations in the one-way model 
(e.g. Hadamards, vr/S gates, and controUed-Z gates), 
and compose the operations for these unitaries to find 
an algorithm for the composite unitary U . We may then 
transform the measurement algorithm so that all of the 
entangling operations are performed first. 

Is it possible to develop measurement algorithms with- 
out reference to the circuit model? One proposal 
reduces the problem of implementing a unitary in the 
one-way model to a problem of expressing the complex 
coefficients of the unitary operator to be implemented in 
terms of sums of roots of unity, which define an entangle- 
ment graph through their ratios. Doing this requires that 
one determine the order in which the measurements are 
to be made. This may be done by making use of a flow 
in the sense introduced in which is a property of just 
the entanglement graph and the vertex sets I,0 C V{G). 

In this article, I describe how to efficiently determine 



whether a graph G (together with input/output vertices 
/ and O) has a fiow in the sense of 0] , for the special case 
/| = \0\ corresponding to algorithms performing unitary 
transformations (as opposed to general unitary embed- 
dings) . This is done by characterizing flows in terms of a 
families of vertex-disjoint paths, and proving that these 
are unique (when they exist) in the case of |/| = \0\. This 
allows flows to be constructed efficiently, when they exist, 
by reduction to solved problems on directed graphs. 



II. PRELIMINARIES 

In this section, we will review the one-way measure- 
ment model, and the concept of a flow which pertains 
to it. For basic definitions in graph theory, readers may 
refer to Diestel's excellent text [7|. 



Notation and conventions. 

For a graph G, we write V{G) for the set of vertices 
and E{G) for the set of edges of G. Similarly, for a 
directed graph (or digraph) D, we write V{D) for the 
set of vertices and A{D) for the set of directed edges (or 
arcs) of D. If x and y are adjacent, we let xy denote the 
edge between them in a graph, and a; — > y denote an arc 
from X to 2/ in a digraph. When a graph G is clear from 
context, we will write x ^ y for the adjacency relation in 
G. We will use the convention that digraphs may contain 
loops on a single vertex and multiple edges between two 
vertices, but that undirected graphs have neither; and N 
will denote the non-negative integers. 



A. The one-way measurement model 

Computations in the one-way measurement model are 
described by a sequence of primitive operations on a set 
of qubits V. Using the notation of [J|, the permitted 
operations are: 
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• entangling operations E^y^,, which perform a 
controUed-Z operatfon on qubits v,w ^ V ; 

• correction operations consisting of X or Z opera- 
tions on single qubits; 

• measurement operations M", which perform a 
measurement of a single qubit v in an orthonormal 
basis of states in the equator of the Bloch sphere. 

The measurements M" may be described by observables 
of the form Ha — 5 (1 — cos{a)X — sin(a)F); this op- 
erator has eigenvectors \±a} oc |0) ± e*" |1}, with \+a) 
having eigenvalue and |— «) having eigenvalue -1-1. The 
operator M" represents measuring the qubit v using the 
projector Ha, and recording the eigenvalue of the result 
as a bit Sy (referred to as the measurement signal) [T^. 
Later correction or measurement operations can depend 
on the value of Sy, which is referred to as classical feed- 
forward of measurement results. 

When using the measurement-based quantum compu- 
tation to perform "quantum-to-quantum" operations (i.e. 
to transform quantum states), we identify two special 
sets of qubits: the set of input gubits I, which arc not 
operated on by a preparation map, and the set of output 
gubits O, which are not measured. The initial state of 
the qubits in / may be arbitrary, and represents the in- 
put of the measurement algorithm; and the qubits of O 
retain a final quantum state, which represents the output 
of the algorithm. A valid algorithm is any sequence of 
the above operations with the following properties: 

(i) each qubit is prepared at most once and measured 
at most once; 

(a) no operation may depend on a measurement signal 
s„ before the qubit v has been measured; 

{Hi) the first operation performed on a qubit v E P = 
y \ / is a preparation map; 

(iv) the last operation performed on a qubit v £ O'^ = 
V \ O is a measurement. 

Algorithms can be condensed by allowing the measure- 
ment operations to depend on the results of previous 
measurements, and by performing all entanglement op- 
erations towards the beginning of the algorithm. Let s 
be a boolean expression: then, using the equalities 

M^X^ =M-(-i)', (lb) 
M^Z', =M,"+-^ (Ic) 

we may postpone all correction operations until the end 
of the algorithm, and perform all preparation/entangling 
operations at the beginning. This allows us to describe 
measurement algorithms in the usual way with the prepa- 
ration of an entangled resource (an open graph state, de- 
pending on the initial state of the input qubits /), with 
measurement and correction operations performed on it; 
and where the angles of measurements may depend on 



the signals produced by earlier measurements. Note that 
measurements of the form M""'"'^ '' are performed with 
respect to the same basis as M", but with the two ba- 
sis elements |±a) interchanged whenever s = 1 (mod 2). 
Rather than changing the angle of measurement by tt de- 
pending on the value of s, we may add the value of the 
expression s to the measurement result Sy to obtain a 
modified result s^. Equivalently, in the algebraic repre- 
sentation of the measurement algorithm, we may substi- 
tute wherever it occurs in a correction of measurement 
dependency with the expression Si, -I- s: this substitution 
is called signal shifting in (3]. 

Thus, without loss of generality, we may describe 
unitary transformations using measurement-based algo- 
rithms where, if all measurement results are 0, no adap- 
tation of the measurement angles is required, nor are 
corrections on the final output system. This prompts the 
question of whether we can design algorithms in terms of 
the behaviour when all the measurement results are 0. 



B. Flows in the one-way measurement model 

Considering the outward differences between the one- 
way measurement model and the circuit model could lead 
to new techniques for developing quantum algorithms, as 
proposed by [5|. However, an apparent obstacle to di- 
rectly devising algorithms in the one-way model is that 
measurements bases and final corrections on the out- 
put qubits may depend on many measurement signals. 
These details are essential, and raises the questions of 
the order in which measurements are to be performed, 
and what measurement dependencies are required. Un- 
fortunately, this complicates any direct understanding of 
how to perform operations in the one-way measurement 
model, without e.g. translating from the circuit model. 

Given a graph, and a collection of measurements which 
yield a particular operation in the special case where all 
measurement results come out 0, the problem described 
above can be solved if we can determine a measurement 
sequence and signal dependencies from the entanglement 
graph itself, along with the sets of input and output ver- 
tices. In such a sequence of operations, we may treat 
each measurement M" as post-selecting the state of the 
whole system so that v is in the state l+a)', should the 
opposite result occur, the final corrections and the inter- 
mediate changes in measurement angles are equivalent 
to performing a correction immediately after the mea- 
surement on V to bring about the state that would have 
arisen had |-|-q) been the result. We may do this if we 
can infer suitable by-product operations for each mea- 
surement: this can be done if the entanglement graph 
has a flow property introduced in Q. 

Definition 1. A geometry [G, I, O) is an entanglement 
graph G, together with subsets I,0 C V{G) representing 
the sets of input and output vertices of a measurement 
algorithm. A flow on (G, /, O) is an ordered pair (/, ^), 
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FIG. 1: Examples of geometries with flows. Arrows indicate 
the action of a successor function / : O'^ — > P, along undi- 
rected edges. Causal orders =;! for each example are given by 
Hasse diagrams (read from left to right). In the right- most 
example, the two vertices a and b are incomparable, i.e. there 
is no order relation between them. 



FIG. 2: A geometry with no flow. Also shown is a particular 
injection / : O'^ — > I'^, and the coarsest pre-order satisfy- 
ing conditions (|2b|) and — see the discussion on page U 
preceding Definition [T] 



with a function / : O'^ — > P and a partial order \X5^ =4 
on V{G), such that the conditions 

X ^ fix) (2a) 
X 4 fix) (2b) 
y fix) =^ X (2c) 

hold for all vertices a; G and y G ViG). We will refer 
to / as the successor function of the flow, and =:<; as the 
causal order of the flow. 

Figs. [1] and [5] illustrate examples of geometries with 
and without flows. The conditions ((2a|) - pc)) are meant 
to capture a simple set of conditions, independent of the 
angles of the measurements, which are sufflcient to de- 
termine how to adapt measurement angles and perform 
corrections to perform unitary transformations. Specifi- 
cally, it captures when the by-product operations for each 
measurement M" can be considered to consist of a single 
X operation on some qubit f{x), and Z operations on 
each qubit y ^ f{x). 

The above is not a necessary condition for a geome- 
try (G, /, O) to permit unitary evolution independently 
of the measurements performed (see [1] for a general- 
ization), but it is sufficient. The partial order x ^ y 
then represents when such a byproduct operation Cx 
for the measurement on x acts non-trivially on y. If 
we insert the operation Cx after the measurement on 
X (deleting the Pauli operation on x itself) in a mea- 
surement algorithm consisting only of preparation maps, 
entanglers, and measurements, these corrections will be 
absorbed into the measurements on the qubit f{x) and 
each of the qubits y ^ f{x). Performing such substitu- 
tions/absorptions for all of the qubits to be measured, 
the order ^ then describes chains of measurement de- 
pendencies [Tg} . 

Thus, having a flow allows us to infer a sequence of 
measurements, and suitable dependencies for those mea- 
surements, providing a solution to the problem described 



towards the beginning of this section. This makes it eas- 
ier to design quantum algorithms directly in the one-way 
measurement model, by obtaining complete sequencess 
of measurement operations from only partial informa- 
tion. This motivates the problem of efficiently determin- 
ing when a geometry has a flow. 

III. CHARACTERIZING FLOWS IN GRAPH- 
THEORETIC TERMS 

In order to determine whether a geometry (G, /, O) has 
a flow, it is useful to understand the sorts of structures 
which are induced or forbidden in G by the presence of 
a flow. We begin with a restriction of the concept of a 
path cover to geometries: 

Definition 2. Let (G, /, O) be a geometry. A collection 
C of (possibly trivial) directed paths in G is a path cover 
of (G,/,0) if 

(i) each ti G ViG) is contained in exactly one path (i.e. 

the paths cover G and are vertex-disjoint); 
iii) each path in C is either disjoint from /, or intersects 

/ only at its initial point; 
iiii) each path in C intersects O only at its final point. 

In the case |/| = |0|, a path cover of (G, /, O) will just 
be a collection of vertex-disjoint paths from / to O which 
covers all the vertices of G. 

For a flow (/, =^), there is a natural connection between 
the successor function / and path covers for the geometry 
(G, /, O), which we capture in the following Lemma: 

Lemma 3. Let if, ^) be a flow on a geometry (G, /, O) . 
Then there is a path cover Vf of (G, /, O), where x ^ y 
is an arc in some path of Vf if and only if y = fix). 

Proof. Let (/, =5!) be a flow on (G, /, O). Suppose that 



f{x) = f{y) for some x,y G O'^. By condition ((2a|) . we 
have y ^ f{y) = /(x); and by condition (pc|) . we have 
X ^ y. Similarly, we have y ^ x, so x = y. Thus / is an 
injective function. 

Define a digraph P on the vertices of G, and with arcs 
X — > f{x) for a; e O"^. Because / is both a function and 
injective, every vertex in P has maximal out-degree and 
maximal in-degree 1. Thus, P is a collection of vertex- 
disjoint dipaths, dicycles, and closed walks of length 2. 
As well, for every arc {x ^ y) G A(P), we have x ^ y; 
by induction, x ^ z whenever there is a dipath from x to 
z in P. Then if a; and z are such that there are dipaths 
from X to z and from z to x, then x ^ z and z ^ x, in 
which case x — z and the dipaths are trivial. Thus, P is 
acyclic, so P consists entirely of vertex-disjoint dipaths. 

Let Vf be the collection of maximal dipaths in P. 
We show that Vj satisfies each of the criteria of Defi- 
nition [21 

{i) Any vertex v which is neither in dom(/) nor img(/) 
will be isolated in P: then, the trivial path on v 
is an element oi Vf. All other vertices are in ei- 
ther doni(/) or inig (/), and so are contained in a 
non-trivial path oiVf. As these paths are vertex- 
disjoint, each vertex is contained in exactly one 
path. 

{ii) Each vertex in / has in-degree 0, and so may only 

occur at the beginning of any path in Vf. 
{Hi) The vertices in P which have out-degree are pre- 
cisely the output vertices O: therefore one must 
occur at the end of every path, and they may only 
occur at the end of paths inVf. 

Then Vf is a path cover, whose paths contain only arcs 
X —>■ f(x), as required. □ 

It will prove useful to discuss functions / which are 
not necessarily the successor function of a flow (f,^), but 
which nonetheless are related to a path cover in the sense 
of Lemma [21 Thus, we will extend our usage of the term 
successor function to include the following definition: 

Definition 4. Let C be a path cover for a geometry 
{G,I,0). Then the successor function of C is the unique 
/ : O'^ — > I'^ such that y = f{x) if and only if x ^ y is 
an arc in some path of C. If a function / : — > /'^ is 
a successor function of some path-cover of {G,I,0), we 
may call / a successor function of (G, /, O). 

In the case where |/| = \0\, the successor function of a 
geometry (G, /, O) is bijective. This allows us to define 
the additional useful terminology: 

Definition 5. Let C be a path cover for a geometry 
(G,/, O) with |/| = \0\. The predecessor function of C is 
the unique g : I'^ — > O*^ such that g(y) = x if and only 
if X — > y is an arc in some path of C. 

Given that the successor function of a flow for (G, /, O) 



induces a path cover, one might think of also trying to 
obtain a flow from a path cover. There is an obvious 
choice of binary relation which we would like to consider, 
which satisfies conditions (|2bp and ([2c|): 

Definition 6. Let / be a successor function for 
(G, /, O). The natural pre-order [l7j =<! for f is the tran- 
sitive closure on V{G) of the conditions 

X ^ X (3a) 
X =^ fix) (3b) 
y ~ f{x) =^ X 4y (3c) 

for aU x,y e V{G). 

Recall from Section III Bl the description of causal or- 
ders for flows in terms of chains of measurement depen- 
dencies arising from byproduct operations: we have x ^ y 
if there is a sequence of vertices x = zqi -^i, • • ■ , zi = y 
such that the byproduct operator for the measurement 
on Zj acts non-trivially on Zj+i for each ^ j < We 
are then interested in when a natural pre-order is anti- 
symmetric, in which case it provides a well-defined order 
for measurements. 

It is easy to show that the natural pre-order =<; for / is 
a partial order if and only if / is the successor function 
of a flow. If ^ is a partial order, it will be the coarsest 
partial order such that (/, =4) is a flow. However, it is 
easy to construct geometries and successor functions / for 
which the natural pre-order ^ is not a partial order. One 
example is the geometry (G, /, O) illustrated in Fig. [21 on 
page[3l with G equal to the cycle Cq = ao&oai^ia2&2ao> 
/ ~ {ao,ai,a2}, and O — {6o:^ii&2}- For any successor 
function / on this geometry, condition ([3c]) forces either 
Oo ^ ai =^ 02 =<; flo or ap )^ ai )^ 02 )^ Oq to hold. 
Because ag, ai, and 02 are distinct, such a relation =^ 
is not antisymmetric, so it is not a partial order. In 
this case, we have not only a cyclic graph, but a cycle 
of relationships induced by condition ([5c|) . The following 
definitions are aimed to characterize these cyclic relations 
in terms of closed walks. 

Definition 7. Let (G, /, O) be a geometry, and T a 
family of directed paths in G. A walk W = uqUi ■ ■ ■ ui 
is an influencing walk [18] for J- if it is a concatenation 
of zero or more paths (called segments of the influencing 
walk) of the following two types: 

• X y, where this is an arc in some path of J-'; 

• X ^ z ^ y, where x ^ z is an arc in some path of 
T and yz e E{G) is not covered by JT. 

A vicious circuit for is a closed influencing walk for J- 
with at least one segment. 

A non-trivial influencing walk W oi T must start with 
an arc in some path of J-] and that of any two consecutive 
edges of W , at least one is an arc in some path of JT. 
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Then, it is easy to see that the decomposition of W into 
its' segments is unique: the initial segment is of the first 
type if and only if the first two edges are arcs of T, is 
of the second type otherwise. The entire walk can be 
decomposed recursively in this fashion. 

Definition 8. Let (G, /, O) be a geometry. A path 
cover C for (G, /, O) is a causal path cover if C does not 
have any vicious circuits in G. 

The two types of segments which build an influencing 
walk correspond to the flow conditions (|2bp and (Pc)) : 
influencing walks again represent chains of dependencies 
induced by byproduct operations. Specifically: 

Lemma 9. Let C be a path cover for (G, /, O) with suc- 
cessor function f, and let =4 be the natural pre-order of 
f . Then x ^ y if and only if there is an influencing walk 
for C from x to y. 

Proof. To show that x ^ y \i there is an influencing 
walk from x to y, we proceed by induction on the num- 
ber of segments of the influencing walk. If the number 
of segments of the influencing walk is zero, then x = y, 
in which case x =4 y. Otherwise, suppose the proposi- 
tion holds for all influencing walks for C with fewer than 
n segments for some n € N, and that there is an in- 
fluencing walk W — xui ■ ■ ■ ugy from x to y (for some 
vertex-sequence (ujY^^^) which has n segments. 

• If the final segment of W is ugy, then xui ■ ■ ■ ui is 
an influencing walk of n — 1 segments, so x ^ ug. 
Because we also have y — f{ui), from the definition 
of the natural pre-order we have x ^ y. 

• If the final segment of W is ui-iU£y, then 
xui ■ ■ -Ui-i is an influencing walk of n — 1 seg- 
ments, so X =4 ui^i. Because we also have y ~ 
ug = f{ui-i), from the definition of the natural 
pre-order we have x ^ y. 

Conversely: from the definition of =^ as a transitive 
closure, ifx ^ y for some x,y £ V(G), there is a sequence 
of vertices (uj) for some ^ e N such that mq = x, 
Um = y, and either Uj+i — f{uj) or Wj+i ~ f{uj) holds 
for each ^ j < ^. Then, define the paths 



U]f{Uj)Uj+i , 



if Uj+i = f{uj) 
if Uj+i - f{uj) 



(4) 



for each ^ j < £: the walk gqcti ■ ■ ■ cji obtained from 
concatenating these paths is a walk from x to y, and in 
particular an influencing walk. □ 

This equivalence allows us to characterize flows in 
terms of paths and circuits in the graph: 

Tiieorem 10. Let {G,I,0) be a geometry with path 
cover C, f be the successor function ofC, and =^ be the 



natural pre-order for f . Then the following are equiva- 
lent: 

(i) C has no vicious circuits; 
(ii) ^ is a partial order; 
(Hi) (/, ^) is a flow. 

In particular, a geometry has a flow iff it has a causal 
path cover. 

Proof. By construction, (/, ^) fails to be a flow if and 
only if ^ is not a partial order (i.e. if and only if it is not 



anti-symmetric). Thus {ii 



{Hi] 



If =<; is not anti-symmetric, then there are distinct 
x,y G V{G) such that x ^ y and y ^ x: by Lemma IH 
there is then an influencing walk W with at least one 
segment from x to y, and also an influencing walk W' 
with at least one segment from y to x. Then WW is 
an influencing walk with at least two segments from x to 
itself, and is therefore a vicious circuit for C; then C is 
not a causal path cover. Conversely, if C is not a causal 
path cover, then there is a vicious circuit xuiU2 ■ • • ug-ix 
for C: if U2 = f{vLi), then xui and U1U2 ■ ■ - x are both 
influencing walks, in which case x ^ ui ^ x] otherwise, 
XU1U2 and U2 - ■ ■ ui^ix are both influencing walks for C, 
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in which case x ^ U2 ^ x. Thus (i) 



(n) 



□ 



Characterizing flows in terms of causal path covers al- 
lows us to shift the emphasis from the constructibility of 
a causal order =:<; to the absence of vicious circuits. By 
using successor and predecessor functions, we may show 
that requiring vicious circuits to be absent for a path 
cover yields a strong uniqueness result: 

Theorem 11. Let (G, /, O) be a geometry such that 
\I\ = \0\. If {G,I,0) has a causal path cover C, then 
C is also the only maximum collection of vertex- disjoint 
dipaths from I to O. 

Proof. Suppose that C is a path cover for (G, /, O) with 
successor function / : O*^ — > J'^, and suppose there is a 
maximum-size collection J- of vertex-disjoint / - O di- 
paths which differs from C. Let S C V{G) be the set of 
vertices not covered by IF: because \T\ = \C\ = \I\ = |0|, 
we have S* n / = and 5' fl O = 0, in which case T 
is a path cover for the geometry (G \ S*, /, O). Then, let 
g' : {I'^ \ S) — * {C^ \ S) be the predecessor function of 
T as a. path cover of (G \ 5, /, O). 

Because C and J- differ, there must exist a vertex x G 
such that x f{x) is not an arc in some path of J-. 
Note also that for v e dom(/), f{v) ^ dom(/) holds only 
if f{v) e O \ I C dom(g'); that is, f{v) e dom(/) U 
dom(g'). Then, deflne a vertex sequence {uj)j^j^ in G by 
setting Mo = x, ui = f{x), and 

f f{uj) , Uj eS or Uj / f{uj-i) 
Uj+i = < (5) 
I g'{uj) , Uj ^ S and Uj = 

for all j ^ I. Fig. [3] illustrates this construction. 
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FIG. 3: An influencing walk for a path cover C (solid arrows) 
induced by another maximum collection T of vertex-disjoint 
paths from 7 to O (hollow arrows). The shaded area is a 
subset of the set S C V{G) not covered by T. 



Clearly Uj 



J-j+i 



for all j G N. We also have uq 



Ul 



Because Wm' is an 

-1 Urn. or Urn. ^ 



an arc in some path of C, and for any j 1 such that 
Uj — > Uj+i is not an arc of C, it follows that Uj ^ f(uj-i), 
in which case we have Uj^i — f{uj), which implies Uj 
Mj+i is an arc in some path of C. Then for any iV e N, 
the walk Wn = uqUi ■ ■ ■ un is an influencing walk in G. 

Because G is a finite graph, the Pigeon Hole Principle 
implies that there must be integers to, m' G N with to < 

77l', Um = Um', and Um-l = Mm'-l 

influencing walk, at least one of u. 
Um+i is an arc in some path of C. In the former case, 
the closed walk Um-iUm ■ ■ • Um'-i is an influencing walk, 
and thus a vicious circuit; otherwise, UmUm+i • • ■Um' is 
an influencing walk, and thus a vicious circuit. In either 
case, there exists a vicious circuit for C, in which case C 
is not a causal path cover. 

Thus, if C is a causal path cover, there can be no such 
vertex sequence [uj)^^^ as defined above, and so there 
can be no maximum family of vertex-disjoint I - O paths 
T which differs from C. □ 

Note that for |/| = |0|, because a causal path cover of 
(G, /, O) is unique if it exists, and the successor function 
of any flow will also be the successor function of a causal 
path cover, there is at most one successor function / 
which yields a flow for (G, /, O). Because the natural pre- 
order ^ for / is coarser than any other valid causal order 
for /, it too is unique. Then, a "mimumum-depth" flow 
for a geometry (G, /, O) is unique in the case |/| = \0\. 



IV. FINDING FLOWS EFFICIENTLY WHEN 

|/| = |o| 

Theorem 1111 allows us to reduce the problem of finding 
a flow for (G, /,0) when |/| = \0\ to finding a maxi- 
mal collection C of vertex-disjoint I - O paths, and then 
determining whether or not C has vicious cycles. Both 
steps can be expressed in terms of solved problems in di- 
rected graphs, and both can be solved in time 0{km), 
where fc = |/| = \0\, and m — \E{G)\. An upper bound 
on the number of edges that a geometry may have if it 



has a flow allows us to further bound this by O(fc^n), 
where n = \V{G)\. In this section, I give an outline for 
the solution of these results to show that a flow can be 
found efficiently when |/| = |0|. [H 



A. Finding a path cover for {G, I, O) 

Finding a path cover for (G, /, O) can be reduced to an 
instance of network flows. A network is a directed graph 
N with a designated source vertex r and sink vertex s, 
and a capacity function c : A{N) — > N representing the 
maximum rate at which some substance can pass through 
each arc. An integral r - s network flow is a function 
(j) : A{N) — > N such that (/)(a) < c(a) for all a e A{N), 
and where the "net flow" into a vertex x £ V{N), 



$(x) = 



eA{N) 



eA(jv) 



(6) 



is zero for x ^ {r, s}. The value of the network flow cj) is 
$(r). 

We may start the reduction to network flows by aug- 
menting the entanglement graph G to a graph G', adding 
a vertex r which is adjacent to every vertex of /, and a 
vertex s which is adjacent to every vertex of O. Any col- 
lection of vertex-disjoint I-O paths then corresponds nat- 
urally to a collection of "internally disjoint" paths from 
r to s of the same size. By a construction presented in 
Section 8.3 of fTl|, we can then efficiently construct from 
G' a network N with source r and sink s, such that every 
integral r - s network flow </> can be used to construct a 
collection of $(r) internally disjoint paths from r to s in 
G'. It then suffices to find a maximum integral network 
flow for N. This is a well-studied problem: in the case 
where all edges have capacity 1, the Ford-Fulkerson algo- 
rithm (see e.g. [l^. Section 26.2) runs in time O(fc'm'), 
where k' ^ |/| = \0\ is the value of the maximum net- 
work flow, and where m' = 0{\E{G)\) is the number of 
arcs in the network N. 

Having found a maximum-size collection J- of vertex- 
disjoint paths, we may determine if is a path cover 
simply by verifying that it covers all vertices: this may 
be done in time 0(|T^(G)|). If is not a path cover, 
(G, /, O) has no flow by Theorems [TO] and [IH 



B. Determining a causal order 

To determine whether or not a path cover C (with suc- 
cessor function /) for (G, /, O) has vicious circuits, we 
may create the digraph Jf whose vertices are those of G, 
and where (x — *■ y) € A(J7/) iff there is an influencing 
walk for C of at most one segment from x to y. Then, 
C has vicious circuits iff J7/ contains a directed cycle. 
Tarjan's algorithm (see e.g. [l^. Section 3.1) is a simple 
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algorithm for determining the strongly connected compo- 
nents of a directed graph D: the equivalence classes of 
vertices which are mutually reachable by directed walks. 
In any circuit of J7/ , all of the vertices are mutually reach- 
able; then, we can use Tarjan's algorithm on J7/ to de- 
termine whether C is a causal path cover. If J7/ contains 
two mutually reachable vertices, (G, /, O) has no flows 
by Theorems [TU] and [TlJ 

Because the natural pre-order ^ for / is characterized 
by influencing walks for C, we have x ^ y iff there is a 
directed path from a; to y in J7/. Then, the problem of 
computing =<! is equivalent to the problem of computing 
the transitive closure of J7/: the directed graph Tf in 
which there is an arc from a; to y iff there is a non-trivial 
directed walk in J7/ from x to y. The transitive closure 
can also be computed by a modification of Tarjan's al- 
gorithm: then, =^ can computed at the same time as we 
are determining whether J7/ contains directed cycles (i.e. 
whether or not ^ is anti-symmetric). 

Each path of C is totally ordered by the pre-order =:^: 
then, we can represent the relation ^ eflticiently through 
a chain decomposition — for each x G V{G), we store the 
minimal element j/p in each path P £ C such that x ^ yp. 
From Theorem 3.11 of [l3l|, we can compute =<; in time 
O(fcm), where fc = |C| = |/| = \0\ and m = \E{G)\. As 
remarked in the previous paragraph, we may determine 
whether ^ is a partial order at the same time: if it is, 
(/, is a flow for (G, /, O); otherwise, (G, /, O) has no 
flows. 



C. Eliminating geometries with too many edges 

For n — \V{G)\ and k — \0\ (but without requiring 
that / and O have the same number of elements) , an ex- 
tremal result shows that any geometry (G, /, O) which 
has a flow has at most kn — C^^^) edges. Thus, if G has 
more than this number of edges, it cannot have a flow. 
We can use this as a preliminary test for any geometry 
when deciding if it has a flow, aborting if G has too many 
edges: if a geometry passes this test, the subroutines for 
finding the successor function / and the partial order =^ 



described above run in time O(fc^n). 



V. SUMMARY AND OPEN QUESTIONS 

Flows are a tool for analyzing the underlying geome- 
try of measurement algorithms, which may make it fea- 
sible to develop algorithms in the one-way model with- 
out direct reference to the circuit model. We have seen 
how they can be characterized and efficiently found using 
tools of graph theory, in the special case where the input 
and output systems have the same number of qubits. 

One direction in which this work could be general- 
ized is by considering the generalization of flows (called 
"gflows") presented by Browne, Kashefi, Mhalla, and 
Perdrix [8] describe a which accommodate more complex 
byproduct operations for measurements. An efficient al- 
gorithm for finding such gflows would be a substantial 
advance in the line of investigation of flnding appropri- 
ate byproduct operations. 

Another question is whether similar work can be done 
in "classical-to-classical" measurement-based quantum 
computing, wherein all qubits are prepared and mea- 
sured. Is it possible to find sequences for measure- 
ment, and appropriate measurement dependencies, so 
that complete algorithms for universal quantum compu- 
tation (including the final measurements, and without 
initial states other than |-|-)®") can be obtained in the 
one-way measurement model only from partial informa- 
tion (the entanglement graph and the observables to be 
measured)? 
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